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ABSTRACT 

We present detailed spatially resolved measurements of the thermodynamic properties of the 
X-ray emitting gas in the inner regions of the five nearest, X-ray and optically brightest, and 
most X-ray morphologically relaxed giant elliptical galaxies known. Beyond the innermost 
region at r > 1 kpc, and out to r ~ 6 kpc, the density, pressure, entropy, and cooling time 
distributions for the X-ray emitting gas follow remarkably similar, simple, power-law like 
distributions. Notably, the entropy profiles follow a form K oc r^, with an index a = 0.92-1.07. 
The cumulative hot X-ray emitting gas mass profiles and the gas-mass to stellar-light ratios of 
all five galaxies are also similar. Overall the observed similarity of the thermodynamic profiles 
in this radial range argues that, in these systems, relativistic jets heat the gas at a similar 
rate averaged over time scales longer than the cooling time ^cooi ^ 10^ yr. These jets are 
powered by accretion from the hot gas, or material entrained within it, onto the central super- 
massive black hole. This jet heating creates an energy balance where heating and cooling are 
in equilibrium, keeping the hot galactic atmospheres in a 'steady- state'. Within r < 1 kpc, 
this similarity breaks down: the observed entropy profiles show well resolved flattening and 
the values difl'er from system to system substantially. The accretion rate onto the black hole 
and the AGN activity, heating the interstellar medium, must therefore vary significantly on 
time scales shorter than ^cooi = 10^-10^ yr. 
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1 INTRODUCTION 

X-ray studies with Chandra and XMM-Newton have shown that 
relativistic jets, produced by accreting supermassive black holes 
(SMBH) at the centres of galaxies, groups and clusters, interact 
strongly with their environments, driving shocks (e.g.|Forman et al. 
2005 2007 , Nulsen et al.|20 05', "Simionescu et al. 2009b , Million| 
et al. 2010b ) and inflating bubbles of relativistic plasma in the sur- 
rounding X-ray emitting gas (e.g. |Churazov et al.| 2000[ [ 



[ |Fabian| 



eFar'200?, ^2006^, "Birzan et al.'"2004; 'Dunn et al. '"2005 '/Dunn &' 
Fabia n 2006 2008 , Forman et al. 2005 2007 , Raff^erty et al. 2006 , 
McNamara & Nulsen|2007p . These bubbles rise buoyantly and can 



entrain and uplift large quantities of low entropy gas from the inner- 
most regions of their host galaxies ( [Simionescu et al.|2008||2009a| 
[Werner et ari[20T0l [20TI) . All of this activity is believed to take 
place in a tight feedback loop, where the hot interstellar medium 
(ISM), which is observable in the X-ray band, cools and accretes 
onto the central SMBH (referred to as active galactic nucleus or 
AGN), leading to the formation of jets which heat the surrounding 
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gas and drive the above phenomena (for a review see [McNamara[ 
[& Nulsen[[2007[ [Gitti et al.|2012| ). This heating acts to lower the 
accretion rate, reducing the feedback, until the accretion eventually 
builds up again. AGN feedback appears to play an important role in 
regulating the 'cooling cores' of galaxy clusters and the formation 
and evolution of massive galaxies (e.g. [Peterson & Fabian[[2006[ 
[Croton et al.|2006|[DeLucia & Blaizot|2007[[Sijacki et al.|2007| ). 
Although much has been learned about AGN feedback, many ques- 
tions regarding the remarkable balance between heating and cool- 
ing remain. 

Nearby giant ellipticals are in many respects low redshift prox- 
ies for more distant and more luminous cluster cooling cores, al- 
lowing us to study AGN feedback 'in closeup'. Another key advan- 
tage is that the binding energy per particle is lower in galaxies: a 
given amount of AGN heating therefore aff'ects the gas in galaxies 
more obviously than in clusters. By observing a sample of nearby, 
X-ray bright galaxies with a high degree of completeness, and with 
data of sufficient quality to map the detailed thermodynamic prop- 
erties, we can hope to gain important insight into the physics of 
AGN feedback, its duty cycle, and its role in galaxy formation. 

One of the most useful thermodynamic quantities in studying 
the impact of cooling and feedback in galaxies, groups, and clusters 
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Table 1. Summary of the Chandra observations. Columns list the distances of the galaxies, the angular scales per kpc at the corresponding distances, the 
B-band optical and bolometric X-ray luminosities jO' Sullivan et al. 2001), radio lumin osities at 1.4 GHz jCondon et al.|1998||2002) , the line-of-sight Galactic 
absorbing hydrogen column densities, Nu, in their directions jKalberla et al.|2005) , observation dates and identifiers, detectors, and exposure times after 
cleaning. 



vj cilci A y 




Scale 








iVH 




Obs. ID 








(Mpc) 


(arcsec kpc~^) 


(101%) 
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NGC4472 (M49) 


16.7* 


12.4 


8.7 


2.96 


1.20 


1.53 


2000 Jun 12 


321 


ACIS-S3 


19.1 








5.8 








2010 Feb 27 


11274 


ACIS-S3 


39.7 


NGC4649 (M60) 


16.5* 


12.5 




2.05 


0.13 


2.04 


2007 Jan 30 


8182 


ACIS-S3 


49.5 
















2007 Feb 01 


8507 


ACIS-S3 


17.5 


NGC1399 


20.9* 


9.9 


4.4 


5.68 


1.52 


1.50 


2000 Jan 18 


319 


ACIS-S3 


49.9 
















2003 May 26 


4172 


ACIS-I 


36.7 
















2008 Jun 08 


9530 


ACIS-S3 


59.3 


NGC1407 


28.8"^ 


7.2 


7.4 


1.95 


1.23 


5.42 


2000 Aug 16 


791 


ACIS-S3 


32.5 


NGC4261 


31.6'f^ 


6.5 


5.1 


1.63 


329.52 


1.75 


2008 Feb 12 


9569 


ACIS-S3 


100.9 
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iBlakeslee et al.A2009] 



of galaxies is the gas entropy which, by rewriting the adiabatic in- 
dex in terms of X-ray observables, can be defined as ^ = kTn~^^^, 
where kT is the temperature and is the electron density of the 
X-ray emitting gas. Because entropy can only be changed by gains 
and losses of energy, it traces the thermal history of the gas. The ra- 
dial profiles of entropy in statistical samples of galaxy clusters and 
groups have been analyzed b yjPratt et al.| { [2006l|2010| ), [Cavagnolo| 
[eTaLl ( [20091 ) and |Sun et aLl ( |2009| ). However, a detailed study of 
the thermodynamic properties of giant ellipticals, reaching down to 
sub-kpc scales, has not been performed. 

Here we present detailed, spatially resolved measurements of 
the key thermodynamic properties (density, temperature, pressure, 
entropy, cooling time, and gas mass) of the X-ray emitting gas 
in the innermost regions of the five nearest. X-ray and optically 
brightest, and most morphologically relaxed giant elliptical galax- 
ies known. All of these galaxies have deep Chandra observations, 
as well as exquisite optical and multi- wavelength data. 

In Sect. [2] we describe our sample, the reduction of the X-ray 
data, and the spectral analysis. In Sect.|3] we present the measure- 
ments of both the 2D distributions and the deprojected radial pro- 
files of the key thermodynamic properties. We show, for the first 
time, that the azimuthally averaged gas density, pressure, cooling 
time, and entropy profiles of relaxed giant elliptical galaxies are in 
the radial range 1 ^ r < 6 kpc remarkably similar. Based on these 
results, in Sect. [4] we present some new insights into the likely ther- 
mal history of the X-ray emitting atmospheres of the galaxies. In 
Sect.|5] we summarize our conclusions. 



2 SAMPLE SELECTION, DATA REDUCTION, AND 
ANALYSIS 

2.1 Sample selection and its properties 

Our target list is drawn from the parent sample of jPunn et al.| 
( |2010| ), who identified the optically-brightest (Bj ^13.5 mag) and 
X-ray brightest (Fx > lO'^^ erg cm'^ s"^ in the 0.1-2.4 keV band) 
giant elliptical galaxies within a distance J ^ 100 Mpc and with 
declination dec ^ -45 degrees. Motivated by a goal to resolve the 
innermost thermodynamic structure in detail, from this list we se- 
lect only those systems with distance d < 35 Mpc. From the result- 
ing sample of thirteen galaxies, we select those systems that appear 



the most morphologically relaxed at X-ray wavelengths (i.e. sym- 
metric and undisturbed in their Chandra images). The highly sym- 
metric morphologies of these galaxies ensure that their deprojected 
thermodynamic properties can be determined as robustly as pos- 
sible. Five galaxies meet these criteria. Their X-ray, near-infrared, 
and radio images (except for NGC 1407, for which we do not have 
spatially resolved radio data) are shown in Figs.[TJ|5]and their prop- 
erties are summarized in Table [T] 

As seen in the X-ray images (the reduction of the X-ray data 
is described in Sect. [2^, NGC 4649 and NGC 1407 are highly 
symmetric. Although NGC 1399, NGC 4472, and NGC 4261 do 
display slight departures from spherical symmetry, they are still 
among the most X-ray morphologically relaxed giant elliptical 
galaxies known. The faintest galaxy in our sample is NGC 4261, 
with a bolometric X-ray luminosity of Lx = 1.63 X 10^^ erg s~^ 
and flux of 1.36 x 10~^^ erg s~^ cm~^. Our relatively high lumi- 
nosity limit restricts the contamination of the detected X-ray emis- 
sion by an unresolved population of X-ray binaries (see Sect. |2.2.2] 
for details). The relatively high X-ray fluxes of our sources, com- 
bined with their proximity ensure that we can measure the thermo- 
dynamic properties down to sub-kpc radii. 

The processed near-infrared K-band images obtained with the 
Two Micron All-Sky Survey (2MASS) shown in the upper central 
panels of Figs. [TJ|5] were dow nloaded from the NASA/IPAC In- 
frared Science Archiv ( Jarrett et al.||2003t . They show that the 
stellar light and by implication the star-mass distribution in all of 
the galaxies is highly symmetric with no indications for disturbance 
or star-formation. The galaxies do not, in general, show detectable 
dust lanes or H i emission and their stellar populations are old and 
red ( |Serra & Oosterloo 12010] ). NGC 4261 does show a nuclear disk 
of gas and dust in the inn ermost r ~ 100 pc (see |Jaflre et al.|199^ 
JaflTe & McNamar a'1997) containing (5.4 + 1.8)xlO'^ M© of molec- 
ular and atomic hydrogen ( Ferrarese et al.||1996| ); however, cool 
gas or dust has not been detected elsewhere in NGC 4261 ( jSerra &| 
|Oosterloo|2010[ [Combes et al.|2007| ). All five galaxies are essen- 
tially 'red and dead' . 

All five galaxies have central, active radio jets, but their ob- 
served radio powers and morphologies span a wide range, as in- 
dicated by the radio luminosities listed in Table [T] and by the ra- 

^ http://irsa.ipac.caltech.edu/ 
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Figure 1. Top panels: Chandra X-ray image of NGC 4472 in the 0.5-2 keV band (left panel); K-band near-infrared image from the Two Micron All-Sky 
Survey (2MASS; central panel); VLA 1.4 GHz radio map obtained in the C-configuration from jPunn et al.|{2010) with the X-ray contours over-plotted (right 
panel). The images are shown in on a logarithmic scale. Bottom panels: 2D map of pressure (in units of keV cm~^ x (20^) ^^^^ panel), entropy (in units 

of keV cm^ x (20!^)^^^' central panel), and temperature (in units of keV; right panel). The maps were obtained by fitting each region independently with a 
single temperature thermal model, yielding Icr fractional uncertainties of ~3-5 per cent in temperature, entropy, and pressure. Point sources were excluded 
from the spectral analysis and therefore appear as black circles on the maps. 



dio maps in Figs. [T}{4] (see also |Dunn et"aLl[20T0l ). NGC 4261 is 
the most powerful radio source, with jets that puncture through the 
X-ray emitting galactic atmosphere to form giant lobes in the sur- 
rounding intra-group medium (IGM). NGC 1399 has ~10 kpc long 
radio jets with lobes that also appear to be depositing energy away 
from the centre of the galaxy. The radio lobes of NGC 4472, and 
of the 10 times less radio luminous NGC 4649, are seen at smaller 
radii. For NGC 1407, we do not presently have high quality radio 
data, but its radio luminosity is similar to those of NGC 4472 and 
NGC 1399 (see Table[TJ, clearly indicating the presence of jets. 

2.2 Chandra data reduction and analysis 

2.2.1 Data reduction 

Our analysis of the Chandra data follows the data reduction proce- 
dures described in Million et al.| ( |2010a ) and Mill ion et al.| ( |2010b| ). 
The standard level- 1 event lists were reprocessed using the CIAO 
(version 4.3) software package, including the latest gain maps and 
calibration products. Bad pixels were removed and standard grade 
selections applied. The data were cleaned to remove periods of 
anomalously high background. The net exposure times after clean- 
ing are summarized in Table IT] Background images and spectra 



were extracted from the blank- sky fields available from the Chan- 
dra X-ray Center. These were cleaned in an identical way to the 
source observations, reprojected to the same coordinate system, and 
normalized by the ratio of the observed to blank- sky count rates in 
the 9.5-12 keV band. 

Background subtracted images were created in 6 narrow en- 
ergy bands, spanning 0.5-2.0 keV. These were flat fielded with re- 
spect to the median energy for each image and then co-added to 
create the broad band X-ray images shown in Figs. [TJ|5] 

Identification of point sources was performed using the CIAO 
task WAVDETECT. These were excised from the images in Figs. 
[TJ|5] and substituted by local backgrounds. Point sources were 
also excluded from all regions used for spectral analysis. Separate 
photon-weighted response matrices and eflfective area files were 
constructed for each spectral region. 



2.2.2 Spectral analysis 

Our spectral analysis was separated into two parts. In the first part, 
2D maps of projected thermodynamic quantities were made. The 
individual regions for the 2D spectral mapping were determined us- 
ing the Contour Binning algorithm ( [Sanders |2006| ), which groups 
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Figure 2. As Fig.[T] but for NGC 4649. The VLA 1.4 GHz radio map was obtained in the B -configuration jPunn etal.|2010) . 



neighboring pixels of similar surface brightness until a desired 
signal-to-noise threshold is met. In order to have small enough 
regions to resolve substructure, yet still have enough counts to 
achieve better than 5 per cent accuracy in the temperature mea- 
surements, we adopted a signal-to-noise ratio of 18 (~320 counts 
per region). We modeled the sp ectra extracted from each spatial re- 
gion with the SPE^I^ package |Kaastra et al.|l996| SPEX is faster, 
uses an updated version of the MEKAL plasma model with respect 
to XSPEC, though provides consistent results). Because the spec- 
tral band above 2 keV is dominated by background in low tem- 
perature galaxies, the spectral fitting was performed in the 0.6-2.0 
keV band. The spectrum for each region was fitted with a model 
consisting of an absorbed single-phase plasma in collisional ion- 
ization equilibrium, with the temperature and spectral normaliza- 
tion (emission measure) as free parameters. The line-of-sight ab- 
sorption column densities, Nu, were fixed to the values determined 
by the Leiden/ Argentine/Bonn radio survey of Hi ( [Kalberla et al.| 
|2505l see Table[TJ. 

After the point source removal, the contribution from an un- 
resolved population of low mass X-ray binaries (LMXB) to the 
X-ray emission is insignificant. We performed fits, where we in- 
cluded power-law emission components with photon indices F = 
1.56 - 1.80, which were shown to describe well the population of 
low-luminosity LMXB (Irwin et al. 2003 , Kim & Fabbiano 2004). 
Such additional power-law emission is not detected significantly 
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in the galaxies and including it in the spectral fitting does not sig- 
nificantly aflfect the measured thermodynamic properties of the hot 
ISM (we note that this conclusion does not change if we extend the 
fitted band to 7 keV). 

The largest systematic uncertainty in the measured densities 
stems from the metallicity measurements. When fitting the spectra 
of thermal plasmas with kT in the range ~0. 5-1.0 keV, the emission 
measure and metallicity typically anti-correlate. In regions with 
complicated temperature structure the metallicity is often biased 
low and consequently the emission measure is biased high (e. g. 
|Buote||2000[ [Werner et aL]|2008| ). Moreover, for spectra with rel- 
atively low numbers of counts, the statistical uncertainties on the 
metallicity are high and multi-temperature model fits are unfea- 
sible. Therefore, our 2D maps of thermodynamic properties were 
produced with the metallicity fixed to 0.5 Solar (throughout the 
paper, metal abundances are given with respect to the Solar val- 
ues by [ Grevesse & SauvaT||1998| ). If the metallicity is underesti- 
mated/overestimated by a factor of 2 in the model, the density will 
be overestimated/underestimated by a factor of ~ 1.35. 

In the second stage of our analysis, we measured azimuthally- 
averaged, deprojected radial profiles of thermodynamic quantities. 
We extracted spectra from concentric annular regions at least 1 arc- 
sec wide, with a signal-to-noise ratio of at least 18 (~320 counts). 
The deprojected profiles of thermodynamic properties were ob- 
tained with the XSPEC (version 12.5 |Arnaud| 1996) spectral fitting 
package, using the projct model. The combined set of azimuthally 
averaged spectra extracted from concentric annuli was modeled in 
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the 0.6-2.0 keV band simultaneously to determine the deprojected 
electron density («e) and temperature (^T^e) profiles. The emission 
from each spherical shell was modeled with a photoelectrically ab- 
sorbed ( |Balucinska-Church & McCammon||1992] ) APEC thermal 
plasma model ( [Smith et al.|2001 using AtomDB v2.0.1). Depro- 
jected densities were determined in 14-27 shells, and temperatures 
in 7-10 regions, depending on the galaxy. For each galaxy we fit- 
ted two metallicity values, one for the shells outside the radius of 
~ 1 kpc and another value for the shells in the central regions of the 
galaxies. 

For all spectral fitting, we employ the extended C-statistics 
available in XSPEC and SPEX. All errors are quoted at the 68 per 
cent confidence level for one interesting parameter (AC =1). 



3 THERMODYNAMIC MEASUREMENTS 
3.1 2D distribution of thermodynamic properties 



The lower panels of Figs. [TJj5] show 2D maps of projected electron 

pressures (n^kTo), entropies {kT^lnl), and temperatures (^T^e). The 
projected electron densities were determined assuming a fixed line 
of sight of column depth / = 20 kpc. 

The 2D distributions of the pressures are approximately spher- 
ically symmetric in the cores of all five galaxies, consistent with 
the X-ray emitting gas being in approximate hydrostatic equilib- 
rium. The most obvious departure from spherical symmetry is seen 
in NGC 4472 in the areas spatially coincident with the radio lobes 



at radii r ~ 4 kpc, where the projected thermal pressure drops by 
10-25 per cent. 

The spatial distribution of the projected entropy appears spher- 
ically symmetric in 4/5 systems, the exception being NGC 1399 
where the entropy distribution beyond r ~ 1.5 kpc is elongated, 
with lower entropy gas extended along the radio jets. Overall, how- 
ever, in all cases, the lowest entropy gas resides in the innermost 
cores of the galaxies, as expected for convective stability. The re- 
gions of relatively high entropy at r > 6 kpc, to the north of 
NGC 4472 and NGC 1407, and to the east of NGC 1399, indicate 
that these galaxies are moving through the ambient medium. 

The observed temperature distributions of the galaxies appear 
less symmetric. In NGC 4472 we see cooler plasma extended along 
the eastern radio lobe, and in NGC 1399 the low temperature ISM is 
distributed along the radio jets. The lowest projected temperatures 
of the galactic atmospheres are in the range 0.6-0.75 keV. Such low 
ISM temperatures are not, however, always seen in the densest, cen- 
tral regions. Three out of five galaxies (NGC 4649, NGC 1407, and 
NGC 1399) show a significant temperature increase in their cores. 
The clearest central temperature increase is seen in NGC 4649, ar- 
guably the most morphologically relaxed galaxy in our sample. 



3.2 Azimuthally-averaged, deprojected thermodynamic 
properties 

The relatively undisturbed X-ray morphologies and high degree of 
azimuthal symmetry observed in the distributions of projected ther- 
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Figure 4. As Fig.[T] but for NGC 4261. 



modynamic quantities confirms that the five systems in this study 
are among the most dynamically relaxed giant elliptical galaxies in 
the local Universe. All systematic uncertainties associated with de- 
projection analyses, carried out under the assumption of spherical 
symmetry and of a single-phase nature at each radius, are therefore 
minimized for these systems. We have determined the azimuthally- 
averaged deprojected radial profiles of electron density (Hq), tem- 
perature (^re), and metallicity (Z) from the Chandra data. Using 
these quantities we also determine the deprojected radial profiles 

2 

for the cumulative gas mass, entropy {K - kT^jnl ), electron pres- 
sure (Pe = n^kTo), and cooling time. We define the cooling time as 
the gas enthalpy divided by the energy radiated per unit volume of 
the plasma: 



^cool ~ 



nQniA{T) 



(1) 



where the ion number density rii = 0.92nQ, and A(T) is the cool- 
ing function for Solar metallicity tabulated by |Schure et al. ([2009]) 



Cooling functions based on older plasma codes f Sutherland & Do-| 
ipita 1993 ) predict a 9 per cent lower cooling rate for 0.5 keV 
plasma. 

Our main result from the deprojection analysis is that beyond 
the central r ~ 1 kpc, and out to the radii r < 6 kpc, the cumulative 
gas mass, density, pressure, entropy, and cooling time distributions 
for the X-ray emitting gas follow remarkably similar radial profiles 
(see Fig.|6}. The observed cumulative gas mass within r = 8 kpc 
is ~ 5.8 X 10^ 



cent. The entropy profiles follow a simple, power-law form K oc 
with an index a = 0.92-1.07. The weighted system-to-system 
dispersions in entropy at r = 1 kpc and at r = 5 kpc are only 14 per 
cent and 17 per cent, respectively. The density profiles follow the 
form Hq oc with an index = 1.1-1.46 and with a system-to- 
system dispersion of 10 per cent at r = 5 kpc. The pressure profiles 
have an approximate power-law form of oc with an index 
7 = 0.93-1.43 and with a system-to-system dispersion of 13 per 
cent at r = 5 kpc. The thermodynamic profiles of NGC 4261 can 
be well described with power-law relations all the way down to the 
smallest radii measured. 

In the absence of heating, cooling in the central regions of 
these galaxies would be strong: within the radius r < 4 kpc, the 
cooling times in all five galaxies are ^cooi < 1 Gyr. Within r < 1 kpc 
of the center, the cooling times are ^cooi < 10^ yr, and less than a 
few 10^ yr within 0.5 kpc. 

All five galaxies show a flattening of their central entropy pro- 
files within r < 1 kpc. In contrast to the remarkable similarity 
in the thermodynamic profiles at larger radii (1 ^ ^ ^ 6 kpc), 
within this innermost r < 1 kpc region the entropies of the galax- 
ies also show a significant system-to-system scatter of 40 per cent 
at r = 0.2 kpc. The system with the least central flattening of 
the entropy distribution is NGC 4261, where the entropy drops to 
^0 = 0-88 + 0.02 keV cm^ at 0.23 kpc. The largest core entropy of 
TTo = 4.1 + 0.5 keV cm^ at r = 0.27 kpc is observed in NGC 1407. 

The temperature profiles of the galaxies exhibit significant 
system-to-system scatter. Interestingly the temperature profiles for 
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Figure 5. As Fig.[T] but for NGC 1407. 



NGC 4649, NGC 1399, and NGC 1407 show clear temperature 
increase in the central regions, which may be associated with re- 
peated shock activity increasing the entropy of the gas. In contrast, 
the radial temperature distributions in the centres of NGC 4472 and 
NGC 4261 are relatively flat. 

The best fit metallicities of the hot ISM within r ~ 1 kpc are 
Z ~ 0.5 Solar, except for NGC 1399 which has a core metallicity 
of Z ~ 0.85 Solar. Beyond this innermost region (1 < r < 6 kpc) 
the metallicities of NGC 4649 and NGC 4261 are Z ~ 0.3 So- 
lar; NGC 1407 has Z ~ 0.75 Solar, and the metal abundances of 
NGC 1399 and NGC 4472 are Z ~ 1 Solar. The systematic uncer- 
tainties on the metallicity measurements of < 1 keV plasma 
with CCD type detectors are, however, unfortunately still signifi- 
cant. 

To verify that the spectroscopically measured temperatures, 
and other thermodynamic quantities, determined in azimuthally av- 
eraged annular regions are approximately unbiased with respect to 
the true average values at a given radius, we have compared their re- 
sults with the average values determined from all bins at the same 
radii in the 2D maps. We find no significant diff'erences between 
the values found by fitting spectra extracted from circular annuli 
and the azimuthally averaged values determined from 2D maps at 
any radius. 

The obtained deprojected density and temperature profiles are 
broadly consistent with the previously published results in the lit- 
erature (e.g. [Humphrey et al.|2006| [20091 piurazov al.|2010| ). 



4 DISCUSSION 

Beyond the innermost core, at r > 1 kpc and out to at least 
r ~ 6 kpc, the X-ray emitting hot gas mass profiles and ther- 
modynamic profiles exhibit remarkable self- similarity. The den- 
sity, entropy, pressure, and cooling time distributions follow sim- 
ple, power-law forms. The total gravitating mass at these radii is 
also similar in all five galaxies, and is dominated by their stellar 
populations ( Humph rey et al.|[2006[ [Das et aLl[2010| ). The mass 
fraction of hot X-ray emitting gas in this radial range is small, 
of the order of /gas ^ 2 x 10~^. Because the hot ISM represents 
a minor and essentially dynamically insignificant fraction of the 
total mass of the galaxies, the observed similarity of the hot gas 
mass and its thermodynamic profiles is surprising and reflects a 
deep and long term stability of energy input from the AGN and 
from the stellar populations. Interestingly, the ratio of the cumu- 
lative gas mass within r < 8 kpc to the total B-band luminosity 
for the galaxies in our sample spans a relatively small range of 
Mgas/LB = 7.1 - 11.3 X 10-3 Mo/Lo. 

4.1 Universal thermodynamic profiles 

The similarity of the thermodynamic properties of the gaseous at- 
mospheres of our galaxy sample has important implications for the 
evolution of giant ellipticals. The cooling time of the hot gas in 
the central regions is short (see left panel of Fig. |6]) and in the ab- 
sence of heating this gas will cool and form stars (e.g. [Peterson &| 
|Fabian|2006| ). The X-ray spectra of these galaxies, however, lack 
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the expected strong O vi i lines associated with g as cooUng below 
r ~ 5 X 10^ K (see e.g. [Sanders & Fabian|20rT| . Moreover, their 
stellar populations are old, showing no evidence for recent star- 
formation, and no significant central atomic or molecular reservoirs 
of cooled gas are detected ( Serra & Oosterloo|2010| ). All of this in- 
dicates that the gas has been kept from cooling. 

Type la supernovae are expected to explode in g iant ellipti- 
cal galaxies at a rate of 0.166 per 100 yr per 10^^ Lb© ( [Cappellaro 
|et al. 1999 ), releasing ~ 10^^ ergs of energy per supernova. In prin- 
ciple, this can provide a heating rate of 1.6 x 10"^^ erg s~^ within 
r < 6 kpc, which is of the same order of magnitude as the X- 
ray luminosity within the same region. However, to approximately 
balance the radiative cooling, the energy released by supernovae 
would have to couple to the X-ray emitting gas with a 100 per cent 
efficiency, which seems unlikely. The dominant source of heating 
is instead likely to be the central AGN. 

As indicated by the observed radio morphologies, the jets and 
lobes driven by the central AGN interact with and mechanically 
heat the ISM out to a radius of several kpc. A complete, flux- 
limited study of 18 nearby, X-ray bright, optically bright giant el- 
liptical galaxies by [Dunn et al.| ( [20T0| shows that 17/18 systems 
exhibit some form of current radio activity associated with the cen- 
tral AGN (see also [Best et al.|2005| [Dunn & Fab ian 20061 [2008) 
|Sun|200 9 ). Such observations imply that active 'radio-mode' AGN 
feedback is not a rare and sporadic phenomenon, but rather rep- 
resents the default state for large elliptical galaxies. [Allen et aT] 
( |2006| ) demonstrate a tight correlation between the Bondi accretion 
rates of hot gas and the observed jet power in such ellipticals. This 
argues that AGN feedback is likely controlled by the accretion of 
hot gas, or multiphase gas entrained within it. The similarity of the 
observed thermodynamic profiles reported here supports this, and 
argues that relativistic jets, produced by accretion onto the central 
SMBH, heat the gas at a similar and constant rate averaged over 
time scales longer than the cooling time, ^cooi ^ 10^ yr. This jet 
heating creates an energy balance where heating and cooling are in 
approximate equilibrium, keeping the hot galactic atmospheres in 
a near 'steady- state'. Importantly, the observed 'steady- state' also 
indicates that the time averaged accretion rate onto the SMBH is in- 
timately connected to the conditions of the hot X-ray emitting gas 
phase on spatial scales of several kpc. 

The jet-powers estimated from the energies and timescales re- 
quired to inflate the bubbles of relativistic radio-emitting plasma 



are Pjet ~ 8 x 10^^ ergs s'^ in NGC 4472 { [Alien et al. 



Pjet 1.3 X 10^ ^ ergs s"^ in NGC 464 9 and Pjet 



2006| ); 

2.2x10^^ ergs s'^ 

in NGC 139 9 ([Shurkin et al.[[ 2008|); and Pjet > lO^^ ^j-g^ ^-i 
NGC 4261 ((ySullivan et al. 2011). In NGC 4472, NGC 4649, 
and NGC 1399 the energy required to keep the hot ISM within 
r < 6 kpc from cooling is only 2.5-10 per cent of the estimated 
jet-power, and in NGC 4261, which is the system with the most 
powerful radio source in our sample, this energy is less than 1 per 
cent of the power output of the AGN. In NGC 4261 the jets pene- 
trate through the hot galactic corona and form lobes which deposit 
most of their energy into the surrounding intergalactic medium 
( [O'Sullivanetal.|2011| ). 

Although we do not see evidence for shocks currently propa- 
gating through the X-ray atmospheres of the galaxies in our sample, 
repeated AGN induced shocks most likely do contribute to the heat- 
ing of the hot ISM. Detailed studies of the giant elliptical galaxy 
NGC 5813 revealed that the total energy in shocks driven by the 
expanding bubbles can be 1/2-3 times the total internal energy of 
the cavities ( Randall et al.|20H"] ). The total jet power in these sys- 
tems is therefore most likely even higher than the estimates pro- 



vided above. Interestingly, despite the fact that the fraction of the 
estimated 'current' jet-power that is needed in each case to keep 
the hot ISM from cooling varies from galaxy to galaxy by an order 
of magnitude, the observed similarity of the thermodynamic prop- 
erties, including the gas mass and entropy, argues that the energy 
deposited into the ISM within r < 6 kpc must be similar and steady, 
to keep the hot galactic atmospheres in a 'steady-state'. 

4.2 Heating in the innermost r < 1 kpc 



As shown in Sect [3.2| the entropy and cooling time profiles show 
well resolved flattening at radii r < 1 kpc. This flattening, com- 
bined with the short cooling times of ^cooi < 10^ yr, implies that 
the gas at these radii has been heated by AGN activity relatively 
recently. Radiative cooling will reduce the entropy of the central 
gas, and in order to offset such cooling and flatten the profiles a 
heating mechanism is required. Heating by repeated AGN induced 
weak shocks is likely to be the most significant channel for en- 
tropy increase of the hot plasma near to the AGN. Shock heating, 
which increases the entropy by A^S, will off'set the radiative heat 
loss by AQ = TAS (where T is the temperature of the gas). Ex- 
amples of pronounced, ongoing AGN induced shock heating have 
been studied in detail in M 87 (Forman et al. 2005 |2007[[Million[ 
[et al.|2010bi ) and in NGC 5813 ( .Randall etal.2 011). 

The steepest central entropy distribution is seen in NGC 4261, 
where the entropy continues to drop all the way in to the centre of 
the system. The fact that the central entropies at r < 1 kpc, where 
^cooi < 10^ yr, are substantially different from system to system 
indicates that the accretion rate onto the black hole and the result- 
ing jet-power, heating the ISM, varies significantly on time scales 
shorter than a few 10^ yr. 

Given the remarkable similarity of the profiles at r > 1 kpc our 
analysis suggests that the application of hydrostatic mass analyses 
to the regions of those galaxies from 1-6 kpc are likely to be rela- 
tively accurate (e. g. Humphrey et al.|2006{[Churazov et al.|2008| 
2010 ). Hu mphrey et a l. (2008 2009 ) extend their hydrostatic mass 
analyses all the way into the cores of their target galaxies, argu- 
ing that in some systems, most notably in NGC 4649, relative hy- 
drostatic equilibrium applies all the way to small radii allowing to 
infer directly the mass of the central SMBH. The increased system- 
to-system scatter in our observed central gas entropy indicates that 
caution should be employed when extending such analyses to the 
inner 1 kpc, where feedback may induce significant thermodynamic 
variations on timescales as short as a few 10^ yr. 



4.3 The origin of the hot gas 

There are three possible explanations for the origin of the hot gas 
in giant ellipticals: first, it could be leftover baryonic material from 
the process of galaxy formation. In this case, the similarity of the 
gas mass to stellar light ratios in the galaxies would indicate sim- 
ilar star formation and feedback histories. Second, some fraction 
of the X-ray emitting gas could have been accreted from the sur- 
rounding large scale environment. If this process is important, then 
the gas mass to stellar light ratios would depend strongly on the 
large scale environments of the galaxies and would therefore be 
expected to exhibit significant system-to-system variation, which 
is not seen in our sampl^ A third possibility is that stellar mass 

^ The gas mass fractions around the brightest cluster galaxies of cooling 
core clusters are significantly higher than those of the galaxies investigated 
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Figure 6. Left panel: Radial distributions of cooling times, entropies, and cumulative gas masses for our sample of 5 nearby giant elliptical galaxies. Right 
panel: Radial distributions of the deprojected temperatures, electron densities, and electron pressures. 



loss contributes significantly to the X-ray halos (for a review see 
[Mathews & Brighenti||2003| ). In this process, gas from the outer 
layers of evolved stars is mixed and assimilated into the hot phases 
of the galaxies ( |Mathews|1990| ). Hydrodynamic simulations of the 
gas ejected by red giant stars and planetary nebulae have been per- 
formed by Parriott & Bregman| ( |2008j and Bregman & Parriott 
( |2009| ). These simulations predict that about 75 per cent of the 
ejecta produced by red giant stars moving super sonically relative 
to the ambient medium will be shock heated to approximately the 
temperature of the hot ISM. However, the products of stellar mass 
loss evidently do not thermalize universally in all environments. 
The presence of dust and PAHs in the dense cooling cores of some 
galaxy clusters suggests that in those systems much of the ejected 
stellar gas remains cool (jVoit & Donahue 201 1! (Donahue et al.| 
|2011| ). The high ambient pressures in the cores of galaxy clusters 
with massive cooling cores might be promoting more rapid radia- 
tive cooling of the products of stellar mass loss. However, none of 
our galaxies lie at the centers of massive cooling cores. 

While we cannot at present discriminate clearly whether orig- 
inal baryonic matter left over from the process of galaxy formation 
or the products of stellar mass loss dominate the X-ray emitting 
material within 1 < r < 8 kpc, the remarkable similarity of the 
thermodynamic and gas mass profiles, and gas mass to stellar light 
ratios in the systems presented here provides an important, new 
constraint. 



of the five nearest. X-ray and optically brightest, and most morpho- 
logically relaxed giant elliptical galaxies at X-ray wavelengths. Our 
main conclusions may be summarized as follows: 

• Beyond the innermost central region, at r > 1 kpc and out to 
r ~ 6 kpc, the radial profiles of density, pressure, cooling time, and 
entropy of these galaxies follow remarkably similar, simple, power- 
law like distributions. The entropy profiles follow a form K oc 
with an index a - 0.92-1.07, with a system to system dispersion 
at a given radius of only ~15 per cent. The cumulative hot X-ray 
emitting gas mass profiles and the gas-mass to stellar-light ratios of 
all five galaxies are also similar. 

• The observed similarity of the thermodynamic profiles in this 
radial range argues that, in these systems, relativistic jets heat the 
gas at a similar rate averaged over time scales longer than the cool- 
ing time ^cooi ^ 10^ yr- This jet heating creates an energy bal- 
ance where heating and cooling are in equilibrium, keeping the hot 
galactic atmospheres in a 'steady-state'. 

• The entropy profiles show well resolved flattening at radii r < 
1 kpc. In contrast to the remarkable similarity at larger radii (1 < 
r < 6 kpc), the central entropy value diff'ers from system to system, 
with a scatter of 40 per cent at 0.2 kpc. The accretion rate onto the 
black hole and the AGN activity, heating the ISM, must therefore 
vary significantly on time scales shorter than ^cooi = 10^-10^ yr. 

• Our results support the picture in which the jets associated 
with the central AGN are powered by the accretion of hot gas, or 
the material entrained within it. 



5 CONCLUSIONS 

We performed a detailed spatially resolved study of the thermo- 
dynamic properties of the X-ray emitting gas in the inner regions 



in this paper (e. g. [Allen et al.|2008) . Most of the gas in those systems likely 
originates from the surrounding intra-cluster medium. 
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